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O ■ Abstract 

We consider numeration systems where digits are integers and the base is an alge- 
f-H ! braic number /3 such that > 1 and /3 satisfies a polynomial where one coefficient 
^ \ is dominant in a certain sense. For this class of bases /3, we can find an alphabet 
r-| ' of signed-digits on which addition is realizable by a parallel algorithm in constant 
' time. This algorithm is a kind of generalization of the one of Avizienis. We also dis- 
cuss the question of cardinality of the used alphabet, and we are able to modify our 
algorithm in order to work with a smaller alphabet. We then prove that /3 satisfies 
CN ! this dominance condition if and only if it has no conjugate of modulus 1. When the 
base /3 is the Golden Mean, we further refine the construction to obtain a parallel 
00 ' algorithm on the alphabet {—1, 0, 1}. This alphabet cannot be reduced any more. 
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1. Introduction 



A positional numeration system is given by a base and by a set of digits. The 
^ \ base /3 is a real or complex number such that > 1, and the digit set ^ is a finite 
alphabet of real or complex digits. The most studied numeration systems are of 
course the usual ones, where the base is a positive integer. But there have been 
also numerous studies where the base is an irrational real number (the so-called 
P- expansions), a complex number, or a non-integer rational number, etc. Some 
surveys can be found in 12, Chapter 7] and [s]. Chapter 2]. 



In this work, the base (3 is an algebraic number, that is to say, the root of 
a polynomial with integer coefficients, and the digits of A are consecutive integers. 
We consider only finite /3-representations, see Section [2] for definitions. Let us denote 



Fin 



a{P) = { XI U C Z, / finite, Xi G A^ 
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In general, the set Fin^(/3) is not closed under addition. In this paper, we will prove 
that when /3 is a root of a polynomial with a dominant coefficient, the alphabet A 
can be chosen in such a way that Fin^(/3) is a ring, and, moreover, addition can be 
performed by a parallel algorithm with time complexity 0{1). 

Addition of two numbers in the classical 6-ary numeration system, where b is 
an integer ^ 2, has linear time complexity. In order to save time, several solutions 
have been proposed. A popular one is the Avizienis signed-digit representation 
which consists of changing the digit set. Instead of taking digits from the canonical 
alphabet {0, . . . ,b — 1}, they are taken from a symmetric alphabet of the form 
{—a, . . . , 0, . . . , a}, a being an integer such that b/2 < a ^ b — 1 {b has to be ^ 3). 
Such a numeration system is redundant, that is to say, some numbers may have 
several representations. This property allows one to perform addition in constant 
time in parallel, because there is a limited carry propagation. A similar algorithm for 
base 2 has been devised by Chow and Robertson [sj, using alphabet { — 1, 0, 1}. Here 
addition is realized in parallel with a window of size 3. Notice that Cauchy \^ already 
considered the redundant representation in base 10 and alphabet {—5, . . . , 0, . . . , 5}. 

In symbolic dynamics, such functions computable in parallel are called local, 
more precisely p-local, which means that to determine the image of a word by a p- 
local function, it is enough to consider a sliding window of length p of the input. It is 
a general result that a p-local function is computable by an on-line finite automaton 
with delay p — 1, see |6| for instance. 

Amongst the non-standard bases, special attention has been paid to the complex 
ones, which allow to represent any complex number by a single sequence (finite or 
infinite) of natural digits, without separating the real and the imaginary part. For 
instance, it is known that every complex number can be expressed with base + i 
and digit set {0,1}, jlj. 

Parallel algorithms for addition in bases —2, i\/2, 2i and + i have been 
given in 13|. Results on addition in bases —b, iy/b and —l + iin connection with 



automata theory have been presented in [6|. In particular, if b is an integer, \b\ ^3, 
and A = {—a, . . . , 0, . . . a} with a = [■^^'2'^] , addition in base (3 = \fb is computable 
in parallel and is a (g + l)-local function. If |6| ^ 2 is even, set a = \b\/2; then 
addition in base /3 = -^6 on ^ is a {2q + l)-local function. 

The paper is organized as follows: First we suppose that the base /3 is a root of 
a polynomial with integer coefficients such that one of them is greater than twice 
the sum of the moduli of the other ones. We then say that (3 satisfies the strong 
representation of zero property (or, for short, that /3 is SRZ). We give a parallel 
algorithm. Algorithm I, for addition in that case. Theorem 13.21 When /3 is an 
integer ^ 3, it is the same algorithm as the one of Avizienis. 

Section m is devoted to reduction of the working alphabet. We now suppose that 
/3 is a root of a polynomial such that one coefficient is greater than the sum of the 
moduli of the other ones. We then say that /3 satisfies the weak representation of zero 
property (or, for short, that /3 is WRZ). We give a parallel algorithm. Algorithm II, 
where the alphabet is reduced compared to Algorithm I, but there is a fixed number 
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of iterations depending (only) on the weak polynomial satisfied by /3, to be compared 
with Algorithm I which always has just one iteration. 

We then show that in fact all algebraic numbers with no conjugate of modulus 1 
are SRZ (or WRZ), and we give a constructive method to obtain the suitable poly- 
nomial (strong or weak representation of zero) from the minimal polynomial of /3, 
Proposition 15.11 

In the end, we concentrate on the case where /3 is the Golden Mean. Algorithm II 
works with alphabet {—3, . . . , 3}, which is big compared to the minimally redundant 
alphabet. Using ideas similar to the Chow and Robertson algorithm, we obtain a 
parallel algorithm on {—1, 0, 1}, Algorithm III. This algorithm can be applied to the 
Fibonacci numeration system as well. 

2. Preliminaries 

A finite word on the alphabet ^ is a concatenation of a finite number of letters 
of A. The set of words on A is the free monoid A* . The set of bi-infinite words on 
A is denoted A^. 

A finite [5-representation of x with digits in ^ is a finite sequence {xijm^i^ni with 
Xi in A, such that x = Yll=m^iP'' ■ denoted by the word 

•^ji-^n— 1 ■ ■ ■ •^m+l-^m 

with the most significant digit at the left hand-side, as in the decimal numeration 
system. 

The notion of a function computable in parallel comes from computer arithmetic 
(see 10|), where it is defined on finite words, but we give here a definition on bi- 
infinite words. A function ip : — )■ B'^ is said to be computable in parallel if there 
exist two non-negative integers r and t, a positive integer /c, and a function $ from 
A^ to with p = r+t + k, such that if m = {ui)i^i G A^ and v = {vi)i^i G then 
V = (p{u) if, and only if, for every i in Z, Vki+k-i ■■ - Vki = ^{uki+k+t-i ■ ■ ■ Uki-rj^. 
This means that the image oi uhj ip is obtained through a sliding window of length 
p. Such functions are computable in constant time in parallel. 



The notion of a local function comes from symbolic dynamics (see 11|) and is 
often called a sliding block code. It is a function computable in parallel with k = 1. 
The parameter r is called the memory and the parameter t is called the anticipation 
of the function p. The function p is then said to be p-local. 

To be self-contained, we recall the classical algorithms for parallel addition of 
Avizienis [2] , and of Chow and Robertson . In what follows, "for each i in parallel" 
means that "each numbered step is executed in parallel, and the results of the parallel 
computations are shared between the steps" . 



Careful! Indices of Z are decreasing from left to right. 
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Algorithm of Avizienis: Base P = b, b ^ 3 integer, parallel addition on alphabet 
A = {-a, . . . , 0, . . . , a}, 6/2 < a ^ 6 - 1. 

Input: two words Xn - ■ - Xm and Vn - ■ -ym of A* , with m ^ n, x = ^^^m-^*/^* 
Output: a word ■ ■ ■ 2,„ of ^* such that 

n+1 

z = x + y = ^Zif3\ 

i=m 

for each i in parallel do 

0. Zi := Xi + Di 

1. if ^ a then qi := 1, := Zi — b 

if ^ —a then := —1, Tj := + 6 

if —a + l^Zj^a — 1 then := 0, := 

2. 2;^ := + 

Addition realized by the Avizienis algorithm is a 2-local function, with memory 
1 and anticipation 0. Notice that the minimally redundant symmetric alphabet is 
obtained with the value a = \^^~\ ■ 

The Chow and Robertson algorithm works for base 2 and alphabet { — 1,0, 1}. 
We give here a generalization to an even base (3 = b = 2a, A = {—a, . . . , 0, . . . , a}. 

Algorithm of Chow and Robertson: Base P = b = 2a, a ^ 1 integer, parallel 
addition on alphabet A = {—a, . . . ,0, . . . ,a}. 

Input: two words Xn - ■ -Xm and yn - ■ ■ y-m of A* , with m ^ n, x = Y^=m^iP^ ^"^^ 

y = Y.Lmyi(^'- 

Output: a word Zn+i ■ ■ ■ Zm of A* such that 

n+1 

z = x + y = ^Zil3\ 

i=m 

for each i in parallel do 

0. Zi:=Xi + yi 

1. if a + 1 ^ 2j ^ 6 then qi := 1, rj := Zi — b 

if —b^Zi^ —a — 1 then := —1, rj := Zi + b 
if —a + 1^2j^a — 1 then qi := 0, rj := Zi 
\f Zi = a and Zi-i > then qi := 1, := —a 
if 2;j = a and ^ then qi := 0, rj := a 
if Zi = —a and Zj_i < then gj := —1, r^ := a 
if Zi = —a and Zj_i ^ then gj := 0, rj := —a 

2. Zi := qi-i + Ti 
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Addition realized by the Chow and Robertson algorithm is a 3-local function, 
with memory 2 and anticipation 0. 

The main difference between the two algorithms is that the Avizienis algorithm 
is neighbour-free, while the Chow and Robertson algorithm is neighbour-sensitive, 
since the decision taken at position i in Step 1 depends also on the digit at position 
i - 1. 

3. Parallel algorithm 

Now we move to bases /3 that are algebraic numbers. Let us formalize the 
assumption we want /3 to satisfy. 

Definition 3.1. Let /3 be such that > 1. We say that /3 satisfies the strong 
representation of zero property (or, for short, that /3 is SRZ) if there exist integers 
hk, bk-i, . . . ,bi,bo,b_i, . . . ,b_h, for some non-negative integers h and k, such that f3 
is a root of the polynomial 

S{X) = bkX'^ + bk-iX'^-' + --- + b,X + bo + fe-iX-i + ■ ■ ■ + fe_,X-'^ (1) 

and 

bo>2 Yl l^'l • 

i&{-h,...,k}\{0} 

The polynomial S* is a said to be a strong polynomial for /3. 

If (3 satisfies ([1]), then the word b^bk-i ■ ■ -feifeo^-i ■ " " ^-h is a /3-representation 
of zero. From this we can derive a set of rewriting rules, generated by the rule 
bkbk-i ■ ■■bibob-i ■ ■ ■ b^h ^ 0. 

To simplify the notation, we set B = bo and M = ^ \bi \ . The inequality 

ie{-h,...,k}\{o} 

from Definition 13.11 now reads 

B >2M. (2) 
Suppose that f3 is SRZ. We choose the symmetric alphabet 

A = {—a, . . . , 0, . . . , a}, where a = + 

For this fixed alphabet A, we describe a parallel algorithm for addition in base /3. 
Let us denote 

and a'=\^]. (4) 

Then a = a' + cM. The alphabet A' = {—a', . . . , 0, . . . , a'} C A will be referred to 
as the inner alphabet. 



Algorithm I: Parallel addition for base /3 with the strong representation of zero 
property (/3 is SRZ). 



2(B-2M) 



M . 



(3) 



B-l 



2{B-2M) 
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Input: two words Xn - ■ - Xm and Un - ■ -ym of A* , with m ^ n, x = Yll^=m-'^iP^ ^^"^ 
Output: a word ■ ■ ■ z^-h of ^* such that 

n+fc 
i=m— h 

for each i in parallel do 

0. Zi := Xi + y-i 

1. find gj G {— c, . . . , 0, . . . , c} such that Zi — qiB G 

k 

j=-h 



Before proving the correctness of the algorithm, let us stress that the numbers 
c, a, and a' have been defined in such a way that they satisfy 



2a' + 1 ^ B, 



a' + cM ^ a, 



2a-cB ^ a'. 



(5) 



The first two inequalities are clear. To verify the last one, we will use the fact that. 



for any positive real numbers 7, 5, one has \'~f6'] ^ [7] [5] . Putting 7 



B-1 



S = B - 2M, we obtain 



2{B-2M) 



and 



B-1 



2(B-2M) 



[B - 2M) 



B-1 



2{B-2M) 



{B - 2M) = c{B - 2M) . 



This inequality, together with our choice a = a' + cM, already proves the third 
inequality of ([5]). 

Theorem 3.2. Suppose that (3 is SRZ. Then Algorithm I realizes addition in con- 
stant time in parallel in Fin^(/3) with A = {—a, . . . ,0, . . . , a] , where a = + 



B-1 



M. 



2{B-2M) 

Proof. The proof is structured in two parts: 

• The digits of the output belong to A: 

After Line of Algorithm I, Zi belongs to {—2a, 



0, 



2a}. 



Since the inner alphabet A! is formed by 2a' + 1 consecutive integers and 
B ^ 2a' + 1, there exists qi such that Zi — q^B G A! . Since \zi\ ^ 2a and 
2a — cB ^ a', the integer q^ can be found in range {— c, . . . , 0, . . . , c}. 
Finally, in Line 2, the new digit Zi (denoted for the moment by z^'^'^) is 

-1 k 
z^qiB^-J2 ^i-jbj - ■ 

j=-h j=l 
^ ■' 



As \E\ ^ a' and |F| ^ cM, the resulting ^f^'" satisfies \z'^'''"\ ^ a' + cM ^ a. 
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• The output represents the number x + y: 

In order to avoid a tedious description of ranges of summation indices, we will 
consider all coefficients which do not play any role in our consideration to be 
equal to zero. With this convention, fcj = for all i ^ {— /i, . . . , 0, . . . , /c}, and 
therefore we may express the rewriting rule (II]) as ^ = 0. Similarly, Xi 

and yi are set to for all i ^ {m, . . . , n}, and thus x + ?/ = ^ (xj + Also 

the auxiliary coefficients qi are set to for all % ^ {m, . . . , n}. 

After Line 2 of Algorithm I, we obtain 
Since for the last sum we have 

E (E ^^-^^^ /^^ = E E = E E = ' 




the output is correct. 

□ 

Remark 3.3. The correctness of the algorithm stems from the inequalities (ISl), 
which are satisfied by many triplets a, a', c. Our specific choice of a, a', and c 
displayed in ([3]) and (jll) gives the smallest possible value a, and thus the smallest 
cardinality of the alphabet A in our algorithm. Of course, this does not exclude 
existence of other parallel algorithms with smaller alphabets. Let us demonstrate 
the optimality of our choice of a, a', and c. According to (ii) of ([5]), in order to 
minimize a, we have to minimize integers a' and c. The choice of a' in (jlj) is minimal 
with respect to the requirement 2a' + 1 ^ B. By (ii) and (iii) of ([S]), we obtain 

2a' + 1 ^ 2a - 2cM + 1 ^ a' + cB - 2cM + 1 =^ a' ^ cB - 2cM . 

Combining the last inequality with (i) of (E]), and insisting on positivity of c, we 
have 

^ 2a' + 1 ^ 2cfi - 4cM + 1 =^ c ^ 2{b-2M) ' 

which already shows that our choice of c in (jlj) was minimal as well. Let us point 
out that we also used the inequality B > 2M. 

Corollary 3.4. If /3 is SRZ, then addition realized by Algorithm I is a {h+k+l)-local 
function with memory k and anticipation h. Algorithm I is neighbour-free. 
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Example 3.5. Consider the base /3 = 10. It is SRZ for the polynomial — X + 10 = 0, 
where B = 10 and M = 1. In this case 

c = [^] =1, a = I"!] = 5, and a = 6 . 

Therefore, our algorithm provides parallel addition in the decimal numeration sys- 
tem on alphabet A = {—6, . . . , 0, . . . , 6}, and in fact it is precisely the algorithm 
that Avizienis gave in 1961, see |2[. 

More generally, for any integer 6^3, the base P = b is SRZ for the polynomial 
—X + b = 0, and therefore we can apply the same algorithm for addition, namely 
on the alphabet {—a, . . . , 0, . . . , a} with a = \^^~\ ■ 

Example 3.6. Consider the base /3 = 2. For such a base, —X + 2 = is not a strong 
polynomial. Nevertheless, this base satisfies also the polynomial —X'^ + A = 0, which 
already is strong, with B = 4 and M = 1. Now we have 

c=[|]=l, a'=[|]=2, and a = 3. 

So Algorithm I works for base 2 with the alphabet {—3, . . . , 0, . . . , 3}, and is 3-local. 
Remind that the Chow and Robertson algorithm is 3-local as well, but it works with 
smaller alphabet {—1, 0, 1}. 

Example 3.7. Let us consider the base /3 = ^-"y^, the Golden Mean. This base /3 is 
one root of the quadratic equation = X + 1, the second root is (3' = ~ 
Since -|- {P'Y = 7, /3 is a root of the strong polynomial 

S{X) = -X^ + 7-^ 

with B = 7 and M = 2. This implies c = 1, a' = 3, and a = 5. The alphabet we 
work with is .4 = {—5, . . . , 0, . . . , 5}. The inner alphabet is A' = {—3, . . . , 0, . . . , 3}. 

In the following, we denote the signed-digit —b by b (this notation is already 
present in the work of Cauchy 0] ) . Below is shown the performance of Algorithm I 
for addition of two numbers from Fin^(/3), namely 

x = 2/3^ + 5/3^-2/3^ + 5/3^-5/3^ + 3 and y = 5/3S + /3^ + 2/3^-2/3^ + 5/3'^-4/3=^ + 5. 

We have xt ■ ■ Xq = 25255003 and ■ ■ • ?/o = 512254005. 

X ^ 25255003 

y ^ 5122 5 4005 

z ^ 5374 10 9008 

All positions i of z for which Zi does not belong to the inner alphabet {—3, . . . , 0, . . . , 3} 
must be modified by adding or subtracting the rewriting rule 100070001, with the 
digit 7 on position i, as follows: 



8 

















7 


4 


1 n 


q 




n 


s 

O 










r\ 

(J 




1 (J 


(J 


(J 


7 





f\ 




f\ 




1 

























1 











7 











1 























T 











7 











T 























1 











7 











1 























T 











7 











T 



























1 











7 











1 


z 




1 


1 


1 


1 


2 
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Hence x^y = z = (5^^ + 13^^ - (3^ - (3^ + 2(3'^ + 'il3^ + h(5^ -2(5^ + (3^ - (3 + 2- {3-^ + 15-\ 

Example 3.8. Consider the complex base [3 = —1 + i. Since = —4, we can 
use the strong polynomial + 4 with B = 4 and M = 1, z.e., we can perform 
parallel addition on the alphabet A = {—3, . . . , 0, . . . , 3}, and this addition is a 
5-local function. This result appeared in [oj. There is also a parallel algorithm on 
alphabet {—2, ...,2} (a la Chow and Robertson), which gives a 9-local function, 
|6(]. On the minimally redundant alphabet!— 1, 0, 1}, there is a parallel addition 
algorithm provided by Nielsen and Muller [13|. 

Example 3.9. Take for base 13 the rational |. This is an algebraic number, which 
is not an algebraic integer. It is known that any non-negative integer has a finite 
expansion in this base on the alphabet {0, . . . , 6}, see 0]. There is a strong polyno- 
mial S{X) = -2X + 7 with 5 = 7 and M = 2. Thus a' = 3, c = 1, and a = 5. It 
is not difficult to see that any integer has a finite representation on the redundant 
alphabet {—5, . . . , 0, . . . , 5}, and, by Algorithm I, addition is realizable in parallel. 



4. Reduction of the alphabet 

It is very difficult to prove that an alphabet used by some algorithm for parallel 
addition is minimal. It is a lot easier to prove that the alphabet is not minimal; 
which is the case for the alphabet we worked with in Algorithm I. Let us start this 
section with an example. 

Example 4.1. Take for (3 the Golden Mean. This base satisfies the equation 

We say that rewriting rule 10301 H- is positively applied at position i if the number 
3 is added to the actual digit occurring at position i and the number 1 is subtracted 
from the actual digits at positions i — 2 and i + 2. Analogously, we define the 
negatively applied at position i. We shall present a parallel algorithm for addition 
on the alphabet 

^={-3,...,0,...,3} 
Let us explain the algorithm less formally. 
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Input: two finite sequences of digits (xj) and {i/i) of {—3, . . . , 3}, witli x = 
and y = Yl ViP'- 

Output: a finite sequence of digits (zi) of {— 3, . . . , 3} sucli that 



z = x + y = ^Zil3\ 



for each i in parallel do 

0. Zi := Xi + Ui 

1. if 2i e {2,3,4,5,6} 



then apply negatively the rule at i 



if Zi e {-6,-5,-4,- 
if Zi G {-1,0,1} 



3, 



2} then apply positively the rule at i 
then do nothing 



2. if e {2,3,4,5} 



then apply negatively the rule at i 
then apply positively the rule at i 
then do nothing 



if Zi e {-5,-4,-3,- 

if 2, e {-1,0,1} 



2} 



3. if 2, €{2,3,4} 



then apply negatively the rule at i 
then apply positively the rule at i 
then do nothing 



if z, e {-4,-3,-2} 

if z, e {-1,0,1} 



During the course of the algorithm, we change only the actual values of the coeffi- 
cients Zi, but we do not change the sum Y2ziP\ Therefore, for correctness of the 
algorithm, we just have to realize that 

• after Step 0, each Zi belongs to {—6, . . . , 6}, 

• after Step 1, each Zi belongs to {—5, . . . , 5}, 

• after Step 2, each Zi belongs to {—4, . . . , 4}, and, finally, 

• after Step 3, each Zi belongs already to the original alphabet A = {—3, . . . , 3}. 
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Let us illustrate the algorithm on the following example: 



Step 0. 
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Step 2. 
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Step 3. 
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When (3 is the Golden Mean, we have seen how using a weaker representation of 
zero — + 3 — ji = ^ enables to exploit the reduced alphabet {—3, . . . , 3} instead of 
the alphabet {—5, . . . , 5}, which was necessary for applying the strong representation 
of zero —j3^ + 7 — = 0. The idea of this reduction can be generalized to other 
bases (3 as well. 

Definition 4.2. Let /3 be such that > 1. We say that /3 satisfies the weak 
representation of zero property (or, for short, that /3 is WRZ) if there exist integers 
6fc_i, . . . , 6i, ^-1, • • • 5 b^h, for some non-negative integers h and /c, such that (3 
is a root of the polynomial 

W{X) = bkX'' + bk-iX''-^ + --- + biX + bo + 6_iX-^ + ■ ■ • + b_hX-'' (6) 

and 

ie{-h,...,k}\{0} 

The polynomial W is said to be a weak polynomial for /3. 

When P is WRZ, we can describe a parallel algorithm for addition. Let us put 
as above M = X]-t6{-/i fc}\{o} 1^*1' 

^ = {-a, . . . ,0, . . . ,a}, where a = + M . (8) 

Similarly to Algorithm I, the inner alphabet is A' = {—a', . . . , 0, . . . , a'} with a' = 
1"-^^] . The algorithm works in 

s + 1 steps, where s = [-gZjg] • (9) 
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The steps will be indexed by £ = 0, 1, . . . , s. After the •^^'^-step, the digits Zi will 
belong to the alphabet 

= {-2a + e{B - M), . . . ,0, . . . ,2a - i{B - M)} for ^ ^ < s 

and to the alphabet A for i = s. Clearly = {—2a, . . . , —1, 0, 1, ... , 2a}. 

Algorithm II: Parallel addition for base (3 with the weak representation of zero 
property (/3 is WRZ). 

Input: two words x„ ■ ■ - Xm and Un - • -ym of A* , with m ^ n, x = Yll=m.^iP^ ^^"^ 

y = YTi=m.yiP'- 

Output: a word Zn+ks ■ ■ ■ Zm-hs of A* such that 

n+ks 



z = x + y = ^ Zi(3\ 



i=m—hs 



for each i in parallel do 

0. Zi := Xi + y-i 

1. for £ := 1 to s do 

2. if 2;j G A' then := else Qi := sgn Zi 



3. Zi :— Zi Ylij=-h^^-i^. 



Theorem 4.3. Suppose that /3 zs WRZ. Then Algorithm II realizes addition in 
constant time in parallel in Fin_4(/3) with alphabet A = {—a, . . . , 0, . . . , a}, where 
a=\^]+M. 

Proof. During the course of the algorithm, we do not change the value ^ 
(thanks to the fact that we are only applying the weak representation of zero, either 
positively, or negatively). For correctness of the algorithm, we just have to check 
the magnitude of digits at the time when the algorithm stops. 

After the O^'^-step, any digit Zi belongs to the alphabet A^^\ We will prove 

for 1 ^ i < s : if after the {£ — l)*''-step the digit Zi belongs to A^^~^\ then after 
the £*'^-step the new digit Zi (we shall denote it for the moment by z^'^^) belongs 
to A<^'\ 

for i = s : if after the {s — l)*'*-step the digit Zi belongs to A^^~^\ then after the 
s^'^-step the new digit z^'^^ belongs to A. 

This already will confirm the correctness of Algorithm II. 

Let us discuss the value of digits Zi after the f^'^-step according to the value qi 
computed in this step. During the discussion we shall use the inequalities 

a<2a-i{B-M) for i<s, , . 

a^2a- i{B - M) for £ = s, ^ ' 
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which follow from the choice of s by ([9]). Let Zi be in A^^ ^\ There are three possible 
cases: 

• gi = 0: In this case, \zi\ ^ a' = a — M and z^'^^ = Zi — li-j^j- Therefore, 



I ^new I 



/II I I / , ^ r 2a - £(5 - M) if £ < s , 



e for £ < s and 2,"^^" G ^ for ^ = s. 
gi = 1: In this case, a' + 1 = a - M + 1 ^ 2^ ^ 2a - (£ - 1)(5 - M) and 

Therefore, for the upper bound we have 



2a-t{B-M) if £ < s , 
a if £ = s . 



2a + i{B-M) if i<s, 
a if £ = s , 



^ 2a-{£-l){B-M)-B+M = 2a-i{B-M) ^ 
On the other hand, we obtain for the lower bound 
2f ^a-M + l- B- M^-a^ 

i.e., z^^"" G A^^^ for ^ < s and z^^^"" G ^ for £ = s. 
• gj = — 1: Analogous to the previous case. 

In the previous discussion we have also used the inequality a—M^\—B—M ^ —a 
which is a consequence of the definition of a, see ([8]). □ 

Corollary 4.4. zs WRZ, then addition realized by Algorithm II is a {hs+ks+1)- 
local function with memory ks and anticipation hs. Algorithm II is neighbour- free. 

Remark 4.5. If /3 is a root of a strong polynomial, we may use both Algorithms I 
and II with this strong representation of zero. They work with the alphabet A = 

{—a, . . . , 0, . . . , a}, whereby in the first algorithm a = aj = \^-^~\ + 2{b~2M) ^ ' 

while in the second one a = an = l^^l + M. Thus the alphabet of Algorithm II 
is never bigger than the alphabet of Algorithm I using the same SRZ. 

It is easy to see that both the alphabets coincide (a/ = an) if, and only if, 
B ^ AM — 1. If i? and M satisfy this inequality, then the parameter s from 
Algorithm II is equal to 1, and the algorithms are the same. 

If 4M — 1 > i? > 2M, then Algorithm II uses a strictly smaller alphabet. The 
price we have to pay is the number of steps. Unlike Algorithm I, the number of 
steps in Algorithm II depends on values M and B in the weak polynomial. Another 
disadvantage of Algorithm II is the increasing number of positions needed to store 
the result. Comparing the length of outputs in both algorithms, ii d = k + h is the 
degree of the strong or weak polynomial. Algorithm I can enlarge the number of 
used positions by d, while Algorithm II enlarges the number of used positions by ds. 



13 



Remark 4.6. A given base /3 can satisfy more than one strong polynomial, and, 
consequently, we have several versions of Algorithm I. Similarly, one base f3 can 
satisfy several different weak polynomials, and thus several versions of Algorithm II. 

Example 4.7. Let us consider an integer base /3 = 6 in N. 

• lib ^ 3, then —X + 6 = is a strong polynomial for b, with B = b and M = 1, 
thus B ^ 4M — 1. According to the previous remark, the two algorithms (I and 
II) coincide and require the same alphabet A = {—a, . . . , a} with a = • 
In this sense, both algorithms are generalizations of the Avizienis algorithm 
to non-integer bases. 

• For b = 2, the equation —X + 2 = is a weak polynomial and we may use 
it in Algorithm II, which works now with the alphabet A = {—2, —1, 0, 1, 2}. 
This alphabet is bigger than the alphabet in the parallel algorithm given by 
Chow and Robertson, but, on the other hand, it is smaller than the alphabet 
in the algorithm that we have described in Example 13.61 

Remark 4.8. In both our algorithms (I and II), the decision about application of 
the rewriting rule at position i depends only on the actual value of the digit at 
this position. This is the crucial difference from the algorithm described by Chow 
and Robertson for base 2, where the decision depends also on the digit at the right 
neighboring position i — 1. A natural question is whether this kind of strategy may 
be applied to other bases as well. In the last section, we can answer this question 
positively, at least for the Golden Mean. 

5. Polynomials with a dominant coefficient 

Algorithms I and II are applicable only to bases (3 being a root of a strong or 
a weak polynomial, that is to say, a polynomial with integer coefficients, where 
one of these coefficients is dominant, i.e., is greater than twice or once the sum of 
the moduli of the other coefficients. We shall show that for most of the algebraic 
numbers /3 with > 1 such a strong or weak polynomial exists. 

Proposition 5.1. Let a be an algebraic number of degree d with algebraic conjugates 
ai, ^2, . . . , (including a itself). Let us assume that 7^ 1 for all i = 1,2, ... ,d 
and \a\ > 1. Then, for any t^l there exist a polynomial 



Q{X) = aoX™ + aiX""-' + ■■■ + am-iX + G Z[X] 



and an index io G {1, . . . , m} such that 





i£{0,...,m}\{io} 
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Proof. Let G{X) be the minimal polynomial of a, i.e., 

d 

G{X) = - a,) =X'' + g.X"-^ + ■ ■ ■ + g^^^X + e Q[X] . 

i=i 

Let M be the companion matrix of the polynomial G{X): 

I 



M 



-9i 


1 





. 




-92 





1 


. 


. 


-93 








1 . 


. 


-9d^i 








. 


. 1 


-9d 








. 


• 0/ 



^dxd 



It is a well-known fact, which can be easily verified, that the characteristic polyno- 
mial of M satisfies 

det(M-X/) = (-l)'^G(X), 

where / is the unit matrix of corresponding size d x d. In particular, the numbers 
«j are eigenvalues of M. For any n in N, n ^ 2, define 

d 

Gn{X) = lliX - a]) = goin)X^ + gi{n)X''~' + ... + gd^i{n)X + g,{n) . 
Since the matrix M" has eigenvalues a", . . . , a^, we have 



det(M"-X/) = {-lYW{X-a^) = (-1)'^G'„(X) . 

i=i 

As M is a rational matrix, its power is rational as well, and thus det(M"'— X/) is 
a polynomial with rational coefficients. It implies that, for all n in N, the polynomial 
Gn{X) has rational coefficients. Clearly, for all n in N and j in N, j ^ d, we have 



9An) = {-ly E 

{ii,i2,...,ij}£Sj 



n n n 
ii 12 ' ' ' ij ' 



(11) 



where Sj = [A C {1, 2, . . . , rf} : j^A = j} is the set of all subsets of {1, 2, ... , d} 
with cardinality j. 

Without loss of generality, let us assume |ai| ^ \a2\ ^ ... ^ \o:d\, and denote 
jo = max{z : 1 < \ai\}. Our choice of jo guarantees that 



< 1 



aia2 ... a 



Jo 
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for any subset {ii,i2, ■ ■ ■ ,ir} C {1,2, ... ,d} and {ii,i2, . . . ,ir} 7^ {1, 2, . . . , jo}- 
Therefore, for all such choices of {ii, . . . , v} 7^ {1, 2, . . . , jo}, we have 



«1 12 «r 



lim ^^^2 ^ = . 



Since the coefficients gj{n) of the polynomial G„ satisfy ( ITTj) . we can deduce 

lim ^ f for all J = 1, ... , d, j ^ jo , 

a'la^...a]'^ \ {-lY for j = jo . 

We may already claim that there exists uq = nQ{t) in N such that 



or, equivalent ly. 



l^7.oMI>t Yl \9Ano)\, (12) 
ie{i,...,m}\{io} 



JO I je{l,...,m}\{jo} JO I 

The existence of a suitable rio is obvious, as the right hand side tends to and the 
left hand side has the limit 1. 

Let us fix one such uq, and denote by K the least common multiple of the 
denominators of ratios gi{nQ), . . . , Qd-iino), gd{no). Then the polynomial Q{X) = 
KGno{X"'°) and the index io = no jo have the desired properties. □ 

Definition 5.2. Let t ^ 1, and let T{X) = Yl^=o^i-^^ be a polynomial with integer 
coefficients tj in Z. T is called a t-polynomial if there exists io e {0, . . . ,n} such 
that 

\tio I ^ ^ ^ ^ l^i I • 

i6{0,...,n}\{io} 

Remark 5.3. The polynomial Q constructed in the proof of Proposition 15. II clearly 
is a t-polynomial, and, for t = 1 we obtain a weak polynomial, while for t = 2 we 
get a strong polynomial. 



Remark 5.4. Lemma 8 in p| gives a little bit weaker statement in the case of an 
expanding algebraic integer a, i.e., \ai\ > 1 for all conjugates of a. Our proof of the 
previous Proposition 15.11 was strongly inspired by the proof given by S. Akiyama, 
P. Drungilas, and J. Jankauskas in [ij]. 

Example 5.5. The proof of Proposition 15.11 is constructive. The strong represen- 
tations of zero we have used in Examples 13.61 13.71 and 13.81 can be obtained by the 
construction given in the proof with t = 2, namely as follows: 

• For /3 = 2, the minimal polynomial is G{X) = X - 2, then G2{X) =X-2^, 
and the desired Q{X) = G2{X'^) = X^ - A. 
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• For the Golden Mean, the minimal polynomial is G{X) = X"^ — X — 1 = 
(X - f3){X - f3') with (3 = and P' = The minimal n satisfying 

is n = A. Therefore, GiiX) = {X - (3^){X - and the final 

Q{X) = GiiX^) = X^- {(3^ + {P')^)X^ + I3\(3')^ = X8 - 7X^ + 1. 

• For the base /3 = — 1 + the minimal polynomial is G{X) = {X + 1 — i){X + 
1 + = X^ + 2X + 2. Already forn = 2 we have G2(^) = {X + {X - 
(— 1 — i)^) = (X + 2i) (X — 2i) = X^ + 4, and our strong representation of zero 
uses the polynomial Q{X) = Gsl^^) = X^ + 4. 

Proposition 5.6. Let (5 he an algebraic number with a conjugate 7 such that I7I = 1. 
Then /3 cannot satisfy a t-polynomial for any t ^ 1 . 

Proof. Suppose that /3 is a root of a polynomial W{X) = h^X^ + hk-iX^~^ + ■ ■ ■ + 
hiX + 60 + h^iX^^ + ■ ■ ■ + h-hX~^ with integer coefficients hi G Z. Then also the 
conjugate 7 is a root oiW{X): 

hkl^ + hk-i-i^~^ + ■ ■ ■ + 617 + 60 + 6_i7"^ + ■ ■ ■ + h_h-i~^ = 0. 

Since I7I = 1, we obtain 

i£{-h,...,k}\{0} 

□ 

By Proposition 15.11 and Proposition 15.61 we obtain the following result: 

Theorem 5.7. Let P with > 1 he an algebraic number. Then (5 is SRZ (or 
WRZ) if and only if it has no conjugate of modulus 1. 

It is fairly easy to recognize whether an algebraic number does, or does not have 
a conjugate of modulus 1, by looking at its minimal polynomial. First, if the number 
is quadratic, it cannot have any conjugate of modulus 1. Suppose now that a is an 
algebraic number of degree d > 2, with a conjugate a' with modulus \a'\ = 1. Let 
G(X) = X"^ + giX"^-^ + ■■■ + ga^iX + g^ be its minimal polynomial, G{X) in Q[X]. 
Since G{X) is minimal, a' 7^ ±1; i.e., a' is not real. As the minimal polynomial 
has all its coefficients real, the complex conjugate a' = ^ is a root of G as well. In 
general, if the minimal polynomial has two different roots 77 and i, then the minimal 
polynomial satisfies 

GiX)=X'G{^), 

thus it is reciprocal and its degree is even. This is summarized in the following 
remark. 

Remark 5.8. Let /3 with > 1 be an algebraic number of degree d. 

• If d is odd, or 
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• if (i = 2, or 



• if (i ^ 4 is even and the minimal polynomial of /3 is not reciprocal, 
then (3 has no conjugate of modulus 1. 

6. The Golden Mean - all good things come in threes 

In j^, Berstel described an algorithm for parallel addition in base (3 = ^-^-^ on 
an alphabet with cardinality 13. In this section, we give Algorithm 111 for parallel 
addition in this numeration system on the alphabet {—1, 0, 1}. This alphabet cannot 
be further reduced, as proved in [6|. In our algorithm, we use a method similar to 
the method of Chow and Robertson. 

We begin by describing two auxiliary algorithms for elimination of digits. Both 
of them use the weak representation of zero — + 3 — = 0. 

The first auxiliary algorithm removes digits —2: 

Algorithm A: Base (3 = ^^/^ , reduction from alphabet {—2, . . . , 2} to { — 1, . . . , 2}. 

Input: a finite sequence of digits (zj) of {—2, 1, 0, 1, 2}, with z = Zif3\ 
Output: a finite sequence of digits (zi) of {—1, 0, 1, 2}, with z = ^ Zil3\ 

for each i in parallel do 

[ z, = -2 

1. case I Zi = —\ \ then qi 

y Zi = Q and Zi^2 < and Zi_2 < 

else g.j := 

2. Zi := Zi - 3qi + qi+2 + qi-2 



new 
i 



Proof. For the correctness of Algorithm A we have to show that the value z, 
Zi — Sqi + qi+2 + <li-2 belongs to the alphabet {—1, 0, 1, 2} for each i. 
We will use the fact that our prescription for qi satisfies 

Zi and Zi+2 ^ =^ qi = and qi+2 = . (14) 

Let us discuss the value of z'^'^'^ according to the values qi: 

qi = and 

• Zi G {1, 2}: Then, = + qi+2 + q^-2 G {-1, 0, 1, 2}. 

• Zi = 0: Then our algorithm for the value qi supposes that Zi+2 ^ or 
Zi-2 ^ 0. According to ( 1T4|) . we have qi+2 = or gj_2 = 0. Therefore, 

z^"" = qi+2 + q^^2 G {-1,0}. 

qi = —1 and 
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Zi G {-2, -1}: Then, = + 3 + g^+a + g^-2 e {-1, 0, 1, 2}. 

Zi = 0: The rule for the value gj supposes that ^ — 1 and 2;j_2 ^ — 1. 



Therefore, = gi-2 = — 1, which implies zl' 



3 + qi+2 + qi-2 = 1- 



Other combinations, namely (g^ = and Zi G {—2, —1}) or (g^ = —1 and Zi G {1, 2}), 
are impossible. Thereby, we have shown that the digits of the result belong to the 
desired alphabet {—1, 0, 1, 2}. □ 

The second auxiliary algorithm removes digits 2: 

Algorithm B: Base /3 = ^^/^ , reduction from alphabet {—1,0, 1,2} to { — 1,0, 1}. 

Input: a finite sequence of digits (zj) of {—1, 0, 1, 2}, with z = ^ Zi(3\ 
Output: a finite sequence of digits (zi) of { — 1, 0, 1}, with z = Y^ Zi(3\ 

for each i in parallel do 



1. 



case < 



2. 



/ 

Zi 


= 2 










Zi 


= 1 


and {zi+2 ^ 1 


or Zi^2 ^ 


1) 




^z. 


= 


and Zi+2 = Zi^ 


-2 = 2 






Zi 


= 


and Zi+2 = Zi^ 


^2 = 1 and 


2^1+4 ^ 1 


and Zi 


Zi 


= 


and Zi+2 = 2 


and Zi-2 = 


1 and Zi- 


4 > 1 


Zi 


= 


and Zi^2 = 2 


and Zi+2 = 


1 and Zi- 


f4 ^ 1 


qi ■■= 


= 










Zi - 


3gi - 


f gi+2 + gi-2 









> then Qi := 1 



Proof. Before explanation of the correctness of Algorithm B, we point out several 
simple facts about the values g, determined in Algorithm B: 

• Fact 1. Zi ^ 1 and Zi+2 ^ 1 =^ gi = 1 and gi+2 = 1 . 

• Fact 2. Zi ^ and Zi+2 ^ =^ Qi = and gj+2 = . 

• Fact 3. Zi^2 ^ and Zi = 1 and Zi^2 ^ ^> g^ = . 

• Fact 4. Zi = and g^ = 1 =^ qi+2 = 1 and gi_2 = 1 . 
Unlike the previous four facts, the following three ones deserve short proofs: 

• Fact 5. Zi = 1 and g^ = =^ qi+2 = and gj_2 = . 

Proof of Fact 5. The assumption Zi = 1 and gj = implies Zi^2 ^ 0. If 
ZiJ^2 = —1, then we have directly gj+2 = 0. Now suppose that Zj+2 = 0. If 
gi+2 were equal to 1, then, according to Fact 4, the value qi is equal to 1, a 
contradiction. The discussion for gj_2 is analogous. 

• Fact 6. Zi = and g^ = =^ gj+2 = or gj_2 = . 

Proof of Fact 6. When 2;j = 0, then the value qi is zero in these two situations: 
a) One of numbers 2:^+2, Zi_2 is not positive. Then, according to Fact 2, one 
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of the values qi+2, qi-2 is zero. 

b) Both ^ 1 and Zi_2 ^ 1- The assumption qi = forces, without loss 
of generality, that Zi-2 = 1 and Zi-4 ^0. So we can use Fact 3 to deduce 
qi-2 = 0. 

• Fact 7. Zi = 1 and qi = I =^ qi+2 = 1 or gj„2 = 1 • 

Proof of Fact 7. The assumption Zi = 1 and qi = 1 implies 2^+2 ^ 1 or Zi^2 ^ 1- 
According to Fact 1, at least one of numbers gj+2 or qi_2 is 1. 

Now we demonstrate the correctness of Algorithm B, by showing that z"'^"' = 
Zi — 3qi + qi+2 + qi-2 belongs to the alphabet {—1, 0, 1} for each i. 

Let us again discuss the value 2;"^"' according to the values qf. 
qi = and 



Qi 



1 and 



Zi 
Zi 
Zi 

i 

Zi 
Zi 
Zi 



-1: Then, z^^"" 



-l + qi+2 + qi-2 e {-1,0,1}. 
0: Due to Fact 6, we have z^^"^ = gj+2 + qi^2 ^ {0, 1}. 
1: According to Fact 5, we have 2;"'^"' = 1 + qi+2 + qi-2 = 1- 



Then, z^"^"" 



-l + qi+2 + qi-2 e {-1,0,1}. 
Due to Fact 7, we have ^f""' = -2 + qi+2 + qi-2 e {-1, 0}. 
According to Fact 4, we have z^'^^ = — 3 + qi+2 + qi-2 = — 1- 



Other combinations, namely {qi = and Zi = 2) or {qi = 1 and Zi = —1), are 
impossible. Thereby, we have shown that the digits of the result belong to the 
desired alphabet { — 1,0,1}. □ 



Lastly, we can proceed by summarizing the Algorithm III for parallel addition 
in base /3 = ^-y^ and alphabet { — 1, 0, 1}: 



Algorithm III: Base /3 = parallel addition on alphabet A = { — 1, 0, 1}. 

Input: two finite sequences of digits {xi) and {i/i) of { — 1, 0, 1}, with x = J2 ^iP^ 
Output: a finite sequence of digits (zj) of { — 1, 0, 1} such that 

z = x + y = Zi/3\ 

for each i in parallel do 

0. Vi:=Xi + yi 

1. use Algorithm A with input {vi) and denote its output {wi) 

2. use Algorithm B with input {wi) and denote its output {zi) 
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Corollary 6.1. Addition in base the Golden Mean on the alphabet { — 1,0,1} re- 
alized by Algorithm III is a 21-local function with memory 10 and anticipation 10. 
Algorithm III is neighbour-sensitive. 

Proof. It is easily seen that Algorithm A is local with memory 4 and anticipation 4. 
Algorithm B is local with memory 6 and anticipation 6. Therefore, Algorithm III is 
21-local with memory 10 and anticipation 10. □ 

Remark 6.2. All of our algorithms described above are working with symmetric 
alphabets of signed- digits. Such a choice of alphabet is practical, because the al- 
gorithm for addition can be used for subtraction at the same time. However, the 
symmetry of the alphabet is not necessary for parallelism as such. The parallel 
algorithm for addition in base being the Golden Mean given by Berstel uses an 
alphabet of 13 non-negative digits, namely {0, 1 . . . , 9, 10, 11, 12}. In fact, our aux- 
iliary Algorithms A and B demonstrate that, if we do not insist on the symmetry, 
we may be able to further reduce the alphabet. 

Remark 6.3. In Algorithm III, input values X2i and y2i at even positions do not 
influence the output values Z2i+i at odd positions, and similarly for the input in odd 
positions. Therefore, the data at even and odd positions can be stored separately, 
and addition can be performed in parallel on the input (x2i) and {y2i) on one hand, 
and on the input {x2i+i) and (?/2i+i) on the other hand, by an 11-local function 
independently. 

Remark 6.4. The structure of Algorithm III depends heavily on the coefficients of 
the weak polynomial satisfied by the Golden Mean. It is not clear how to generalize 
it to a broader class of bases. However, the case where /3 is a quadratic unit integer 
seems feasible similarly to the Golden Mean. 

Remark 6.5. (Fibonacci numeration system) The Fibonacci numbers are 
defined as follows: 

Fo = 1, Fi = 2, and F„+2 = + F„ forn G N. 

It is well known that any non-negative integer can be expressed as a sum of 
different Fibonacci numbers, 

N = Fi^ + Fi^_^ + . . . + with ik > ik^i > . . . > io ^ . 

If, moreover, we require ij ^ + 2 for any index j, then such an expression is 



unique, see [15[. For example, the number 29 can be represented in the Fibonacci 
numeration system by (29)i;' = 1010000, as 29 = Fg + Ft. If we abandon the require- 
ment of uniqueness, and allow coefficients —1,0,1, then the numeration system is 
redundant, and addition can be performed by a parallel algorithm in constant time 
as well. Let us explain it. We want to add two integers written in the form 

X = ^XiFi and Y = ^yiFi with Xi,yi e {-1,0,1} . 
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We can use Algorithm III with a small modification. Algorithm III is based on the 
weak representation of zero — /S*"*"^ + 3/3* — = 0. The Fibonacci numbers satisfy 
analogous relations: 

• -F^+2 + 3Fi - F,_2 = for alH ^ 2 

• -F3 + 3Fi - Fo = 

• -F2 + 3Fo = 

Therefore, the Algorithm III needs only slight changes, namely at the positions with 
indices i = 0, 1 in Part 1 (Algorithm A) and z = 0, 1, 2, 3 in Part 2 (Algorithm B). 

However, a function computable in parallel is defined by a sliding window of 
length p, which means that the function $ cannot depend on any particular index 
i, but only on the p-tuples of the input alphabet. To get rid of this problem, we use 
the classical trick of extending the alphabet { — 1, 0, 1} by an artificial symbol, say $, 
to indicate the exceptional positions to the function $. Thus we will represent the 
above mentioned number 29 as (29)^ = 1010000$$$ . . .. The domain of the 21-local 
function $ announced in Corollary 16 . 1 1 has to be enlarged to the set { — 1, 0, 1, $}^^, 
and the definition of $ must be correspondingly extended as well. 

Note that, due to the relation — F3 + 3Fi — Fq = 0, the alternative option of 
processing the odd positions separately from the even positions (as described in 
Remark 16.31 for base (3 = ^^/^ ) cannot be used for the Fibonacci numeration 
system. 
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